Spectrum of bound fermion states on vortices in 3 He-B. 
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We study subgap spectra of fermions localized within vortex cores in 3 He-B. We develop an 
analytical treatment of the low-energy states and consider the characteristic properties of fermion 
spectra for different types of vortices. Due to the removed spin degeneracy the spectra of all singly 
quantized vortices consist of two different anomalous branches crossing the Fermi level. For singular 
o and u vortices the anomalous branches are similar to the standard Caroli-de Gennes -Matricon 
ones and intersect the Fermi level at zero angular momentum yet with different slopes corresponding 
to different spin states. On the contrary the spectral branches of nonsingular vortices intersect the 
Fermi level at finite angular momenta which leads to the appearance of a large number of zero 
modes, i.e. energy states at the Fermi level. Considering the v, w and uvw vortices with superfluid 
cores we show that the number of zero modes is proportional to the size of the vortex core. 
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1. Introduction. Since the pioneering work of Caroli, 
de Gennes and Matricon (CdGM) [l[ it is well known that 
quantized vortices in superconductors and Fermi suprflu- 
ids have a non trivial internal electronic structure. It con- 
sists of low energy fermionic excitations localized within 
the vortex cores with characteristic interlevel spacing de- 
fined as Aq/Ef <C An, where A is the energy gap far 
from the vortex line and Ef is the Fermi energy. For 
conventional s-wave superconductors the excitation spec- 
trum of each individual vortex E{Q) of a subgap state 
varies from — Ao to +An as one changes the angular mo- 
mentum Q defined with respect to the vortex axis. 

At small energies \E\ <C A the spectrum is a linear 
function of Q: 



E(Q) ~ -Qu, 



(1) 



where ui « An/(fex£)) Ao is the superconducting gap 
value far from the vortex axis, k± = \Jk F — fef , hp is the 
Fermi momentum, k z is the momentum projection on the 
vortex axis, £ = ftVp/Ao is the coherence length, Vf is 
the Fermi velocity, and Q is half an odd integer. Under 
some exotic conditions Q several vortices can merge and 
then one obtains a multiquantum vortex with a certain 
winding number M. The number of anomalous branches 
per spin projection Q is equal to the vorticity M. For the 
states with an even vorticity all the anomalous branches 
cross the Fermi level at nonzero angular momentum Qf 



E(Q)~-(Q±Q j )A /(k±£) 



(2) 



where j — 1...M/2, Qm/2 ~ k±£. For a vortex with an 
odd winding number there appears a branch crossing the 
Fermi level at zero impact parameter. 

The quantized vortices in 3 He-B have much in com- 
mon with vortices in ordinary s-wave superconductors. 
However in multi-component superfluid system 3 He ax- 
ial symmetry allows the nucleation of additional order 
parameter components inside vortex core. Thus vortices 



in this system are in general nonsingular, i.e. may have 
a superfluid core unlike singular vortices in s-wave su- 
perconductors which always have a normal core. There 
exist five types of vortices with different internal core 
structures in 3 He-B: o, u, v, w and uvw vortices 0, M, @]. 
The o vortex is the most symmetric one, it has no super- 
fluid core and consists of almost pure B-phase without 
inclusions of other phases. Other vortices break some of 
the discrete symmetries existing for the most symmetric 
o vortex. Among them the u vortex is singular while the 
remain v, w and uvw vortices have superfluid cores. 

According to the analysis in the framework of 
Ginzburg-Landau theory [|| @] near the critical tem- 
perature only v vortex is stable. The cores of such vor- 
tices are occupied by an A phase and a ferromagnetic 
phase [1, H, 0| . These additional phases correspond to a 
nonzero total angular momentum projection on the vor- 
tex axis and a zero vorticity in the real space. Therefore 
nuclea of additional phases remain finite at the vortex 
center. Nucleation of ferromagnetic (3 phase inside vor- 
tex cores explains a large spontaneous magnetic moment 
of vortices revealed in the NMR experiments in rotating 
3 He-B Q. The first order phase transition seen in the 
NMR experiments was associated with the change of the 
symmetry of the internal core structure [1, [t| . 

As was shown by Volovik 13, 11 1 in vortices with dis- 
solved core singularity the spectrum of bound fermion 
states can be substantially modified in contrast to or- 
dinary CdGM spectrum of singular vortices. In partic- 
ular the presence of other superfluid phases inside vor- 
tex core leads to the appearance of large number of zero 
modes, i.e. the spectral branches crossing the Fermi 
level. The number of these zero modes can be as high 
as Ef/Aq ~ 3> 1. Thus a minigap in spectrum 
of bound fermions, which is a characteristic feature of 
CdGM spectrum [l| is absent for nonsingular vortices in 
3 Hc-B. Zero modes also exist even for a singular and most 
symmetric o vortex. Although the number of them is 
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much smaller than for nonsingular vortices but it can be 
effectively controlled by external magnetic field II, HI- 
Even in zero magnetic field due to a broken relative spin- 
orbital symmetry in B-phase of 3 He the spin degenera- 



tion of the energy spectrum is removed 12fl . As a re- 
sult the CdGM spectral branches acquire spin depen- 
dent shift which closes the minigap. Localized fcrmions 
which occupy the negative energy states on the spectral 
branches intersecting zero energy level form a one dimen- 
sional Fermi liquid inside vortex core, which can lead to 
the instability of the vortex core structure [l3| ■ 

In this Letter we develop a generalization of CdGM 
theory for the case of vortices in B phase of superfluid 
3 He. We derive a general expression for spectrum of vor- 
tex core quasiparticles in the presence of multiple order 
parameter components inside vortex core and analyze the 
spectra of several particular vortex types. 

The method that we use is based on the approximate 
analytical solution of quasiclassical Andreev equation de- 
scribing the motion of quasiparticles along the trajec- 
tories inside vortex core. Earlier this method was ap- 
plied to study the spectrum of quasiparticles localized 
within the cores of multiquanta vortices 

HQ]. Gener- 
ally the Andreev equation for two component wave func- 
tion ip — (U, V) along the quasiclassical trajectory has 
the form 



~„ dtp ~ 

KT 3 — + T X A R tp 
OS 



f 2 A Im tp = Eip , 



(3) 



where tj. 2,3 are Pauli matrices of Bogolubov-Nambu spin, 
£ is a length scale of the order of coherence length, Ar = 
(A + A+)/2 and iA Im = (A+ - A)/2 are the hermitian 
and anti- hermitian parts of normalized gap operator, E 
is a normalized energy with the normalization energy of 
the order of bulk value of gap function Ao. 

Here we should take into account the spinor structure 
of quasiparticle wave functions which is essential in 3 Hc. 
In this case the coefficients A.Rj rn in Andreev equation 
(|3"|) are 2x2 matrices in spinor space. Then the matrix 
equation l[5]) is a system of 4 scalar equations. If ma- 
trices Ar and Aj m commute [Ar, Ai m ] = the fourth 
order Andreev equation can be reduced to 2 equations 
of the second order. However this can not always be 
the case. To develop a general perturbation theory we 
note we note that if Aj m = the exact solution of 
the Andreev equation ([3]) corresponding to E = can 
be obtained in a spinor basis diagonalizing the matrix 
Ar = diag{ Ari, Ar2)- Then we obtain two degenerate 
solutions "01,2 = (1, ~*)t/i,2(s) corresponding to the zero 
energy, where fj(s) = Aj exp Jq Amds^ and Aj is 

an eigen spinor of matrix Ar. As we will see below in case 
of a single-quantum vortex the functions Ari^{s) have 
asymptotics of different signs A_Ri(+oo)A^i(— 00) < 0. 
We assume that A^+oo) > therefore the solutions 
/i,2(s) decay at s — ±00. Using this localized solution 



as a zero-order approximation for the wave function the 
spectrum can be found within the first order perturba- 
tion theory assuming that \E\ <C I and |A /m (s)| <C 1. 
In general fi t2 (s) are not the eigen spinors of the oper- 
ator Aj m which therefore couples the ipi and ip 2 states. 
Then the standard perturbation theory yields the secular 
equation 
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(4) 



where the matrix elements are <Sii(22) = 
2(A ( 2)|A /m |A (2) ) and S 12 = 2</ 1 |A Jm |/ 2 ). 

In general the accuracy of the first order perturbation 
correction should be determined by the factor 0(A| m ), 
where |A/ m | <C I is a small parameter. However in a 
particular case of Eq.Q the second order correction to 
the zero energy level is exactly zero and therefore the ac- 
curacy of Eq.(j4]) is much better: 0(A 3 m ). To prove this 
result we assume for simplicity that [Ar, A/ m ] = so 
that 5*12 = 0. Then if the eigen function ijj = (U, V) T of 
Eq.Q with Ai m — corresponds to the energy e n the 
other function ip = (— V, U) T corresponds to the energy 
— e„. Therefore it is easy to check that the contribution 
from negative energy levels to the second order pertur- 
bation of the energy level E = exactly compensates the 
contribution from the positive levels. The proof modifica- 
tion to the general case [Ar, Aj m ] ^ is straightforward. 

2. Basic formulas. Our further consideration is 
based on the Bogoulubov- Nambu equation for the quasi- 
particles near the Fermi level. From the beginning we 
assume the system to be homogeneous in z direction 
which coincides with the vortex axis. Then we obtain 
two-dimensional Bogoulubov- Nambu equations with the 
effective Fermi energy E± = EF — h 2 k^/2m and the Fermi 
momentum in xy plane k± — \/k 2 F — k 2 : 



H ip + t x Ar^ - f 2 A Im tp = Etp , 



(5) 



where tj (J — 1,2,3) are Pauli matrices in particle-hole 
space, Ho = ^(p 2 — ti 2 k 2 L )/2m, and p = —ihV. Fur- 
ther we will assume that the gap function and energy are 
normalized to the bulk value of the energy gap Ao . 

Generally the gap function in 3 He-B can be parame- 
terized as follows: A = —ia y (a ■ d), where d is a vec- 
tor in 3D space and cri.2,3 are Pauli matrices in con- 
ventional spin space. Being proportional to the wave 
function of Cooper pairs in isotropic liquid 3 He the gap 
function can be presented as a superposition: A = 
^.^{^"""'''.A^}, where tp is a polar angle 

in xy plane, M is vorticity, A^„ is a Cooper pair wave 
function with definite angular momentum v = —1,0,1 
and spin /i = —1,0,1 projections on the z axis and 
{AB} = (AB + BA)/2 is an anticommutator. 

The order parameter distribution should be axisym- 
metric with the generator of rotation symmetry around 
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z axis [6| Q = L z + S z — MI, where L z and S z are the 
projections of internal angular momentum and spin of 
Cooper pairs onto the z axis and M is a total vorticity. 
Thus for all order parameter components the condition 
(i + v = M should be satisfied. For singly quantized 
vortices M = 1 there can exist five basic components of 
the order parameter. Among them are C\ t —\, C—\ t \ and 
Coo which correspond to the main B phase, Co,i = C a 
and C\fi — Cp which correspond to the additional A and 
(3 phases localized inside vortex core. The additional A 
phase has a zero spin projection (/i = 0) and unit pro- 
jection of orbital momentum [y — 1) on the z axis while 
f3 phase has fj, = 1 and v = 0. The components of ga 
function A^ v are characterized by d vector as follows 
d = X fl (X v ■ q), where q = k/kp and X^ = (x ± iya), 
A = zo • Correspondingly in B phase we have d = q, in A 
phase d = zo(q x + iq y ) and in (3 phase d = q z (x.o + iyo). 
Far from the vortex core at r > only B superfluid 
phase exists so that C\-\ = C-\a = Coo = 1 an d 
Ca,/3 = 0. The vortex type is determined by the be- 



Within the quasiclassical approximation Eq. (|5|) can be 
reduced to 4 equations of the first order along linear tra- 
jectories, i.e. the straight lines along the direction of 
Fermi momentum q = (cos0 p , sin0 p ) (for a detailed re- 
view of this transformation see e.g. Ref . [lij] ) . Each tra- 
jectory is specified by the angle 9 p and the impact pa- 
rameter b — kpZo ■ (q x r). Introducing the coordinate 
along trajectory s = (q- r) we arrive at the quasiclassical 
equation Hip = Eip for the wave function ip(s, 9 P ). The 
quasiclassical hamiltonian is 



d 

H = -i£,q±T 3 — + fiA R - t 2 A It 



(9) 



The impact parameter of quasiclassical trajectories is 
proportional to the projection of angular momentum Q 
of quasiparticles on the z axis: b = —Q/k±. However 
the hamiltonian ([9]) does not commute with the corre- 
sponding operator L z = ~id/d9 p since in general it is 
not conserved in 3 He. Still due to the axial symmetry 
of vortices the total momentum L z + S z is conserved. 



in the quasiclassical hamiltonian (({ 
Let us introduce the new functions 



can be separated. 



U 



'- e " /2 M+U, 



haviour of amplitudes C M „ at smaller distances r ~ f„ Therefore the angular and coordinate variables P and s 
and there exist five types of vortices [6| . 

Vortices of o and u types are singular so that only 
the superfluid components of B phase C\-\, C_i : i and 
Coo are nonzero. These amplitudes are real for the most 
symmetric o vortex and complex ones for u vortex with 
conserved parity P\ — P but broken P 3 = TOJ discrete 
symmetry. Here T is time inversion and is a rotation 
by the angle tt around the axis x perpendicular to the 
vortex axis z. The gap function which describes singular 
B phase vortices can be presented in the form 



(10) 



(11) 



(6) 



where Mo = —ie %Bv a y (a • q). It is easy to check that the 
resulting gap operators (|12ll4ll5p after this transforma- 
tion do not depend on the angle 9 p . For the singular part 
of gap function we obtain 



where d B = (B + q x -iB_q y , B + q y +iB^q x ,C oq z ), B± = 
(Ci _i ± C_i,i)/2. Generally B± = B±(r) are arbi- 
trary complex functions with asymptotics B + (oo) = 1 
and B_(oo) = so that de(oo) = q. 

Nonsingular v, w and uvw vortices have superfluid 
cores with the inclusion of A and f3 phases: 



(s - ib) 



C B + 6A m 



(12) 



where r — Vs 2 + b 2 , Cb = qj^B^ 



q 2 z C 00 and 
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Bl 



q±_ [B- (q z a x - q±a z ) + i{B + - C o)qz^ y ] 



CAqi 



Ap = Cpq z (1 - a z ) 



(7) 



(8) 



(13) 

The expressions for the gap function describing the ad- 
ditional A and (3 phases inside vortex core read 



where q± = k±/kp. The functions Ca,/3 = Ca,p{t) de- 
scribing the spatial distributions of additional A and (3 
phases inside vortex core are finite at r = and van- 
ish outside the core at r > The v and w vortices are 
characterized by real B phase amplitudes. If Ca,(3 are also 
real then we have a v vortex with conserved P2 = PTO x 
symmetry. The case when Re(CA,p) = 0, Im(CA,p) 7^ 
corresponds to w vortex with conserved P3 = TO™ sym- 
metry. The less symmetric uvw vortex with all discrete 
symmetries Pi , P%, P3 broken has complex amplitudes of 
B, A and (3 phases. 



and 



A A = C A {q z q± - iq\v y ) 



A/3 = Cpq z (l — cfg) (q± + q z a x ) 



(14) 



(15) 



Then one can search the solution in the factorized form: 
ip{s,9 p ) = ip(s)exp(i(Q + 1/2)9 P ). Note that from 
Eqs. ljlOllip it follows that the values of azimuthal quan- 
tum number Q should be integer to pertain the unambi- 
guity of initial wave function ip = (U, V) . 

3. Spectrum of vortex core states. At first let 
us consider the quasiparticle spectrum of singular o and 
u vortices when Ca,/3 = 0. In general case Cl_i ^ 
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C_x,i 7^ Co,o the gap function is given by Eqs. (|12ll3|) . 
We will assume the simplifying condition to be fulfilled 
— Cool "C 1 which is justified by Ginzburg- 
Landau calculations @. Then we can take into ac- 
count only the hermitian part of the operator (SAbi)r, — 
(SAbi + 6A B ~ 1 )/2 in Ea.(|12p since the anti hermitian 
part gives the contribution to the energy spectrum of 
the higher order in small parameter 0(|B_|, \B + — Cool). 
In this case the zero order solution of Andreev equation 
(J3j) is spin degenerate ipi^{s) = ipo(s) = (1, —i)rfo{s). 

Diagonalizing the gap function As by spin and using 
the Eq. ([4]) we obtain the energy spectrum in the following 
form 

E(Q,q z ,x) = -u±Q, (16) 

with lo ± = Re(C B /{k^r)) + \F, where F 2 = 
kf [Re 2 (B_/r) Q + q 2 z Im 2 ((B + - C 00 )/r) ] and x = ±1 
corresponds to different spin states. For brevity we have 
denoted (X)o = 2(f \X\f ). For the o vortex with 
Im(B + — Coo) = the difference between w+ and cj_ is 
determined by the asymmetry of amplitudes Ci,_i and 
C-i.i. For the u vortex the condition is less restrictive 
since even in case Ci,_i = C-\ A but C_i : i 7^ Coo we can 
obtain that lu± are different. 

Now we proceed with the analysis of the quasiparticle 
spectra for nonsingular vortices. Here we focus on the 
influence of the additional order parameter components 
and therefore assume the most simple form of singular 
part of gap function (fP2|) with C-l-i — C-x,i — Coo and 
consequently SAbi = 0. At first we consider only the 
influence of A phase and put Cp = 0. In this case we 
have 

(A A ) R = Re(C A )q z q x + Im(C A )q 2 ± a y , (17) 

{A A ) Im = [lm{C A )q ± q z - & y q 2 ± Re(C A )] . (18) 

The matrix coefficients Ajij m in Andreev equation (|3|) 
are diagonalized simultaneously in spinor basis /i^ ~ 
(l,±j) CT so that (Jyfi.2 = ±/i,2- Then Eq.Q yields then 
the following energy spectrum 

E(Q, q z ,x) = -UjQ + ajq ± q z + XJjql, (19) 

where j — 1,2 corresponds to x — 1, — 1, &j — 
2(f 3 \C B /(k±r)\f 3 ) and a hlj = 2Re,Im(f j \C A \f j ). For 
v vortex with Im(C A ) = and w vortex with Re(C A ) — 
it is easy to check that (oj, a, 7)1 = (uj, a, 7)2 in Eq. (fT9|) 
which means following symmetry of v vortex spectrum 
E(Q, q z , x) = —E(—Q, q z , — x). Thus the spectrum con- 
sists of two symmetrical anomalous branches having in 
particular equal slopes as functions of Q at the Fermi 
level E = 0. The spectrum of w vortex is spin divergent 
since jj = in Eq. (ll9p . The situation is more com- 
plicated for the uvw vortex with Re, Im(C A ) 7^ 0. In 



this case (w,a,7)i 7^ (w, 0^,7)2 and the two anomalous 
branches are not symmetrical. Nevertheless there is a 
relation (w, a, 7)i(<? 2 ) = (lo, a, 7)2 ( — 9z) which provides 
spectrum symmetry E{Q,q z ,x) = -E(-Q, -q z , 
corresponding to the CPT invariance of Bogolubov- 
Nambu hamiltonian. 

In Fig. ((TK,b,c) we plot several spectral branches for v, 
w and uvw vortices. We take the model dependencies of 
the amplitudes of B and A components in the following 
form: 

C B {r) =r/^Ttf (20) 

C A (r) = e iK ^/V^Te v , (21) 

where k is an arbitrary phase and are a characteris- 
tic sizes so that at r > ^ the asymptotic behaviour is 
C B (r) -> 1 and C A {r) -> 0. In Fig.JTJI) the two asym- 
metrical anomalous branches sre shown for the uvw vor- 
tex. Note that the CPT invariance is retained with the 
help of anomalous branches corresponding to the oppo- 
site q z values which are shown in Fig.{TJi) with solid and 
dash lines. 

In general case when both A and (3 phases are present 
the spectrum of v vortex still consists of two symmetrical 
branches. When C A ,f) are real the spin structure of zero 
order wave functions is /i^ ~ (q z ,q± =p 1) CT so that the 
hermitian part of gap operator is diagonal in this basis: 

(A.4 + Ap) R = (C A + C p )q z q ± + xqzCp, (22) 

where x ~ il corresponds to fx,2- The anti-hermitian 
part is 

i(A A + Af}) Im = -ia y (q 2 ± C A - q 2 z C ) . (23) 

The spectrum given by Eq. (UJ) differs from Eq. (fTi?]) in the 

absence of j3 pahse: 

E{Q,q z , X ) = -lo+Q + x\/u 2 _Q 2 + \S 12 \ 2 , (24) 

where u>± = (wi±w 2 )/2 and S12 = 2(/ 1 |(A J 4+A ( a)j m |/ 2 ). 
However it is easy to check that spectral branches are 
even functions of q z therefore the spectrum consists of 
two symmetrical anomalous branches as before. 

An analogous procedure for the w vortex with 
Re(C A j3 = 0) yields the spectrum consisting of two spin 
splittcd anomalous branches 

E(Q,q z , X ) = -uQ + 'Y + x\Si2\, (25) 

where uj = ui = lo 2 and 7 = 2(f 1 \(A A + A^)j m |/i>. 
In this case |5'i2 1 (<7z) is even and j(q z ) is odd func- 
tion of q z . Thus the presence of 7 phase removes the 
spin degeneracy of the w vortex spectrum. Two spin 
splitted anomalous branches are not symmetrical, since 
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E(Q,q z ,x) 7^ —E(—Q,q z ,—x) similarly to the case of 
uvw vortex. 

The spectra of bound fermion states of nonsingular 
vortices given by Eqs.fTH I24I25P consist of two anoma- 
lous branches crossing the Fermi level at some points 
Q = Q1.2 7^ 0. Such situation is also realized for the 
spectrum of doubly quantized vortices in ordinary s-wave 
superconductor [see Eq.Q]. An important consequence 
of this fact is an existence of zero modes, i.e. spectral 
branches crossing the Fermi level as one can see in Fig. 
(fT|). Further we discuss zero modes in more detail. 



(a) v vortex: 
Re(C A )*0and lm(C A )= 



(b) w vortex: 
Re(C A )=0 and lm(C A )* 
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FIG. 1: (a)-(c): Several spectral branches corresponding to 
different spin projections \ = 1 (solid lines) and x = — 1 (dash 
lines), (d): Anomalous branches of uvw vortex spectrum. 
Solid and dash lines correspond to the opposite q z values. 



4. Zero modes. As it was shown in Ref. [llj the 

number of zero modes iVo strongly depends on the vortex 
core size £„ . Now with the help of Eq. lfT§)) we will analyze 
this dependence Nq(^ v ) for all three types of nonsingular 
vortices in a model situation when the (3 phase is absent 
Cp = 0. Let us consider the distributions (|2QI2ip of B 
and A order parameter components inside vortex core. 
Such choice of functions Ca,b(^) leads to a simplification 
of the analysis. Indeed in this case we have ctj/ujj — 
q±(k F £ v )s'mn and jj/uJj = q±(k F £, u) cos k. Then the 
right hand side (r.h.s.) of the following equation for zero 
modes does not depend on Q: 

Q 



k F £ v 



= q ± (q z sin k + xq± cos k) . 



(26) 



It is easy to see that for each Q ^ the number of in- 



tersections with Fermi level is two or zero. Thus the 
number of zero modes is Nq — ^{Q m in — Qmax), where 
Qmin/ {k F S,v) and Qmax/(kp£ v ) are minimum and max- 
imum of the r.h.s. of Eq. (|2"6"|) . The additional factor 
of 2 is gained from the summation over the two spin 
states. Then it is obvious that the number of zero modes 
is No — \(k F £ v ), where A ~ 1 is a constant coefficient. 

5. Summary. To summarize we have studied the 
spectra of bound fermion states localized within vortex 
cores for different types of vortices in 3 He-B. In contrast 
to vortices in ordinary s-wave superconductor the spectra 
of singly quantized vortices in 3 Hc-B in general consist 
of two anomalous branches corresponding to the different 
spin structure of quasiparticlc states. This results from 
the removing of spin degeneracy of a standard Caroli- 
de Gennes - Matricon spectrum. The structure of two 
anomalous branches is determined by the vortex type. 

The spectrum of singular o and u vortices given by 
Eq. p^]) is the most similar to the CdGM one |T]). How- 
ever as distinct from the latter it consists of two anoma- 
lous branches with different slopes. As it follows from 
Ginzburg-Landau calculations the asymmetry between 
pairing amplitudes Ci,_i, C-i t i and Coo within vortex 
core is small, therefore the slope difference should also be 
small: |w+ — w_| <C oj±. 

The spectra of nonsingular v, w and uvw vortices con- 
sist of two spin splitted anomalous branches which inter- 
sect the Fermi level at finite values of angular momenta 
Qi,2 7^ 0. In case of v vortex the spectral branches 
are even functions of q z within the same spin subband 
E(Q,q z ) = E(Q,—q z ) which makes the spectrum anal- 
ogous to that of the doubly quantized vortex in s-wave 
supercondictor [see Fig.JTJi)]. 

For w and uvw vortices the spectrum can be a general 
function of q z as it is shown in Figs.([IJ>,c). Note that 
the "skew" of spectral branches of w and uvw vortices 
is produced by the second term in Eq. (|19p which very 
similar to the Doppler shift of the energy a = V s k z which 
would appear due to the superflow V s = ajk F q\ along 
the vortex axis. However there is no real superflow in 
the situation that we consider. Still the order parameter 
symmetry in w and uvw vortices allows the appearance 
of a spontaneous superflow along the vortex axis 0] and 
the effective Doppler shift term in the energy spectrum. 

In contrast to the v vortex the spectra of w and uvw 
vortices (|25|19j) consist of to asymmetrical anomalous 
branches [see Fig.([T])d]. In particular it means that the 
slopes of two anomalous branches dE/dQ can be differ- 
ent at the Fermi level E = 0. As distinct from the case 
of singular o and u vortices the difference in slopes con- 
tains no small parameter and therefore therefore can be 
of significant value. 

Since the anomalous branches in spectra of nonsingu- 
lar vortices intersect the Fermi level at finite angular mo- 
menta Qi.2 =/= there exists a large number of zero modes, 
i.e. the energy states exactly at the Fermi level. We have 
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calculated the number of zero modes for all three types 
of nonsingular vortices assuming a model situation when 
only the additional A phase is present. In a qualitative 
agreement with the results of work 11| the number of 
zero modes was shown to be of the order No ~ kp(,v, 
where is a size of vortex core. 

The significant modification of the spectra of bound 
fermions as compared to the CdGM case should result 
in various ramifications of the vortex dynamics which is 
governed by the kinetics of vortex core quasiparticles [loT ] . 
With the help of the analytical results for the spectra 
obtained in this paper it should be possible to explore 
the dynamics of nonsingular vortices in 3 Hc-B. 
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